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(1967) to evaluate a Coulomb Born matrix element is shown to badly 
misrepresent their total and differential cross sections for the 
Is - 2s transition in atomic hydrogen at low impact energies. 
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PEAKING APPROXIMATIONS IN 
ATOMIC SCATTERING PROBLEMS 

H. Lee Kyle 
M. R. C. McDowell 
NASA, Goddard Space Flight Center 
Greenbelt, Maryland 

hi recent years peaking approximations have been used to evaluate difficult 
integrals in several different models in atomic scattering theory. Vainshtein et 
al. (1963) used one in the Vainshtein approximation. All subsequent calculations 
in the Vainshtein approximation and its modifications (see e.g. Crothers 1967) in- 
clude a peaking approximation. Akerib and Borowitz (1961) and Coleman and 
McDowell (1966) used it in an impulse approximation calculation of direct ex- 
citation by proton impact and Kang and Foiand (1967) used it to evaluate a Coulomb 
Born matrix element in their recent paper. 

In brief a peaking approximation implies the argument that if f(x) is a slowly 
varying function and F(x) is a function with dominant maximum at x = x Q then 

Jf(x)F(x)dx ~ f(x 0 )jF(x)dx . (1) 

Clearly this is a simple version of the method of stationary phase or of steepest 
descents (Morse and Feshbach 1953) but omitting any pretense of mathematical 
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rigor. Use of such an argument will be referred to below as a peaking 
approximation. 

Although this approach has proved useful in calculating radiation induced 
transition probabilities in atoms, see for instance Schiff (1955), its validity in 
connection with most atomic scattering problems is doubtful. The above mentioned 
authors applied it to complicated integrands whose exact behavior is difficult to 
analyze. It is necessary to consider each case individually and until the total 
integrand is exactly evaluated it is impossible to know if the peaking approxima- 
tion is useful in that particular case. 

The Vainshtein cases are being investigated by one of us (H.L.K.) but con- 
clusive results have not yet been obtained. Coleman (1968) demonstrates that, 
for proton impact excitation of atomic hydrogen in the impulse approximation, 
the peaking approximation gives total cross sections which may be as much as 
a factor of 10 in error. In this note we show that in Kang and Foland's theory the 
peaking approximation badly misrepresents the true behavior of the Is - 2s ex- 
citation matrix element for atomic hydrogen. 

Kang and Foland (1967) introduced a new scheme for calculating cross 
sections for the excitation and ionization of atoms and ions by electrons, which 
is essentially an extension of the Coulomb-Born approximation. We do not here 
wish to discuss at length the validity or usefulness of this scheme, but we do 
wish to touch on one point. In dealing with the excitation of neutral hydrogen by 
electrons they represent the incident and scattered electrons by Coulomb waves. 
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These waves do not (asymptotically) properly represent the physical situation. 
Using these Coulomb Born matrix elements they predict finite excitation prob- 
abilities at threshold. The fact is that Coulomb Born matrix elements between 
hydrogenic bound states automatically yield finite excitation cross sections at 
threshold. However one must physically justify the use of such matrix elements 
in the case of neutral atoms. It should be noted that Damburg and Gailitis (1963) 
showed that atomic hydrogen excitation cross sections arising from electron 
impacts are finite at threshold because of coupling effects between degenerate 
quantum states, and not because the scattered electron wave function is Coulombic. 

Our chief concern however is with the peaking approximation used by Kang 
and Foland in evaluating the non-exchange term of their excitation matrix element 
T. f , where 

T i.t = < x . l ' ) (k I )* b (f)lv. b l *.<*>(!ioRa>> • (2) 


Here k Q andkj are the wave vectors of the incident and scattered, a, electron; 
4> b (i) and 4> b (f) are the hydrogenic eigen states of the bound, b, electron, and 
V ab = i/r ab is the interaction potential. The "in" and "out" state, attractive 
Coulomb wave functions of the incident electron are 


*„ <+) U 0 ) 


N 0 e lk °' Ia F(iA 0 , 1. ik 0 r a - ik Q • rj , 
N 1 e lk, ‘ £ ' F(-iAj, 1, - ikj r g - ikj • r a ) 


(3) 


3 



with 


In. I 2 - 2vy / ^k i (l-e 2n/k . 

Kang and Foland show that the direct (non-exchange) excitation amplitude in 
their theory may be put in the form 


T. 


, f 


= N 




i(a-k)-x a 

~k~ 2 


G(k)F(k 0 , r a ) F* (k v rj dr a dk 


(4) 


where G(k) is the generalized oscillator strength, and F(k. , _r ) is one of the 
confluent hypergeometric functions given in (3), and q = k Q -k 1# Since the in- 
tegral with respect to x a increases without limit ask - q , they argue that (4) 
may be replaced by 


T (0) 
i, f 


/ , \ f e i (a-k)-x fl 

= — F(k 0 . r.)F*(k,.r.) dr.dk 


(5) 


This can be put in the form of a Nordseick integral (Nordseick 1954) 


TW = G(q) Nj* N q (l N )x ._ 0 , 


2n Vk 0 /a^ 0 / y + 8 \ ^ \ , i , aS ~ fat \ 

a e \y] \ y I F \ X k Q kj l ' a(y + 8) J 


(6) 
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with 


a - ^ (q 2 +^ 2 ) . fi - k x * q 

y - k 0 * q + iA. kj - a , 8 = k Q k x + k 1 • k Q - /3 , 

0 

and is their version of the peaking approximation to the matrix element, though 
they have used a standard transformation to write the hypergeometric function 
in (6) in an alternative form. Alternative peaking approximations, which differ 
appreciably from (6) may be derived by noting that the slowly varying factor is 
k~ 2 G(k) rather than simply G(k), for example. 

Gailitis (1966) has shown that for s - s transitions (4) may be evaluated 
without approximation, in closed form. We obtain for the Is - 2s function 


T 


Is, 2 s 


4/2 Z 
27 


N o N i 



+ /3 


<321, 


0 <^ 2 J X.= /3 f 


( 7 ) 


with Z being the nuclear charge, and /3 0 = 3Z/2a 0 , and a Q is the Bohr radius. 

We have evaluated the non-exchange contribution to the Is - 2s differential 
cross section in the Kang-Foland model using both (6) and (7). The results for 
the differential cross section 1(8) are shown in Fig. 1, the first Born results 
being given for comparison; (though it is of doubtful value in this energy range.) 
The results are shown as a function of scattering angle 8 - cos -1 (k 0 ‘ kj) at 
energies W = 1.01 and W = 2.0 threshold units (W = 1.0 = 0.75 Ry). Near 
threshold the Kang-Foland method differs markedly from the first Born approxi- 
mation, giving values of 1(8) an order of magnitude larger in the backward di- 


5 


rection. The effect of the peaking approximation is to underestimate by more 
than a factor of 10 at 6 = 0, but to overestimate by 20% in the backward direction. 
At W = 2 the peaking approximation introduces a spurious maximum in I(£) at 
40°, and is generally in error by factors of two or three for 6 < 90°. 

In Table 1 we compare total Is - 2s excitation cross sections obtained from 
(6) and (7) with the Is - 2s - 2p close coupling results [Damburg and Gailitis 
(1963), Omidvar (1964), Burke et al. (1967)] . The Kang-Foland model gives 
results which are approximately a factor of three higher near threshold (W < 1.15) 
than the 3-state close coupling results, while the peaking approximation intro- 
duces an additional error of about a factor of two. At higher energies (W ^ 10) 
both (6) and (7) give total cross-sections in close agreement with the first Born 
approximation. 

We conclude that use of a peaking approximation in the evaluation of matrix 
elements involving Coulomb functions is not in general justified, and may lead to 
order of magnitude errors. All results reported to date in the Vainshtein ap- 
proximation should be treated with suspicion until a detailed analysis of the 
effect of the peaking approximation in that model has been carried out. 
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TABLE CAPTION 


Table I. (Is - 2s) excitation cross sections for atomic hydrogen in units of 

The 3-state close coupling results (c.c.) are from (a) Damburg 
and Gailitis (1963) and (b) Burke, Ormonde and Witaker (1967). The 
Coulomb Born results compare the exact solution of Eq. (7) with its 
approximate solution, Eq. (6), obtained by use of the peaking approxi- 
mation (P.A.). The close coupling calculations of course take account 
of exchange but the other calculations do not. 
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FIGURE CAPTION 


Fig. 1 The effect of the peaking approximation on the differential cross section 
1 ( 0 ) for the Is - 2s transition in atomic hydrogen. Energy is given in 
terms of threshold units (W = 1.0 = 0.75 Ry), (a) W = 1.01, (b) W 
= 2.0. The curves are 1-first Born approximation; 2 -Coulomb-Born, 
Eq. (7); 3 -peaking approximation to the Coulomb-Born, Eq. (6). 
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Table I 


W = k 2/1 

0 ex 

C #C • 

Coulomb Born 

Exact 

P. A. 

1.01 

0.174 (a) 

0.573 

1.13 

1.50 

0.36 (b) 

.396 

-.82 

2.00 

0.25^ 

.301 

.62 
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